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Recently, methanol was identified as a sensitive target system to probe variations of the proton- 
to-electron mass ratio jx [Jansen et al. Phys. Rev. Lett. 106, 100801 (2011)]. The high sensitivity 
of methanol originates from the interplay between overall rotation and hindered internal rotation 
of the molecule - i.e. transitions that convert internal rotation energy into overall rotation energy, 
or vice versa, give rise to an enhancement of the sensitivity coefHcient, Kfj,. As internal rotation is 
a common phenomenon in polyatomic molecules, it is likely that other molecules display similar or 
even larger effects. In this paper we generalize the concepts that form the foundation of the high 
sensitivity in methanol and use this to construct an approximate model which allows to estimate 
the sensitivities of transitions in internal rotor molecules with Cst, symmetry, without performing 
a full calculation of energy levels. We find that a reliable estimate of transition sensitivities can 
be obtained from the three rotational constants (A, B, and C) and three torsional constants (F, 
V3 and p). This model is verified by comparing obtained sensitivities for methanol, acetaldehyde, 
acetamide, methyl formate and acetic acid with a full analysis of the molecular Hamiltonian. From 
the molecules considered, methanol appears to be the most suitable candidate for laboratory and 
cosmological tests searching for a possible variation of fj,. 

PACS numbers: 06.20.Jr, 33.15.-e, 98.80.-k 



I. INTRODUCTION 

Physical theories extending the Standard Model have 
presented scenarios that allow for, or even predict, spatio- 
temporal variations of the constants of nature [1]. Cur- 
rently, a minibcr of laboratory experiments and astro- 
nomical observations are conducted to search for signa- 
tures of such variations. One of the dimensionless con- 
stants that arc hypothesized to vary is the proton-to- 
electron mass ratio, /U = mp/mg. A variation of /U can 
be detected by comparing frequencies of spectral lines 
in molecules as a function of time and/or position. A 
fractional change in fj, will manifest itself as a fractional 
frequency shift. As a measure for the inherent sensitivity 
of a transition, the sensitivity coefficient, K^, is defined 

by 



For pure rotational transitions = — 1, for pure vi- 
brational transitions — while for pure electronic 
transitions if^ = [2] . Transitions between the inversion 
levels of ammonia [3, 4] and hydronium (H30''')[5] have a 
sensitivity of K^^ = —4.2 and K^^ — —2.5, respecively. It 
was shown that the sensitivity of a transition between two 
near-degenerate levels that have a different functional de- 
pendence on fi is enhanced significantly [5 9]. Recently, 
we reported such an enhancement for torsional-rotational 



transitions in methanol [10]. 

Methanol (CH3OH), schematically depicted in Fig. 1, 
consists of an OH group attached to a methyl group. The 
OH and methyl group may rotate with respect to each 
other about the C— O bond. On the left-hand side of 
the figure, the potential energy curve is shown as a func- 
tion of the torsional angle, 7. The interaction between 
the OH and methyl group results in a threefold barrier. 
Tunneling between the three wells, results in a splitting 
of each rotational level into three levels of different tor- 
sional symmetries [11]. Transitions between the different 
torsional levels have a sensitivity coefficient, Kfj, = —2.5. 
As the torsional levels A and E belong to different sym- 
metries, transitions between them are not allowed. It 
was shown in Jansen et al. [10] that transitions convert- 
ing internal rotation energy into overall rotation energy, 
or vice versa, give rise to sensitivity coefficients, /C^, that 
range from -88 to -1-330 in the different isotopologues of 
methanol. 

Hindered internal rotation is a common phenomenon 
found in many polyatomic molecules. Hence, other 
molecules may have similar or larger sensitivities to a 
variation of /x. In this paper we will calculate the sensitiv- 
ities for methanol, acetaldehyde, acetamide, methyl for- 
mate and acetic acid, five relatively small molecules that 
have a group of Csy symmetry that rotates with respect 
to the remainder of the molecule. These five molecules 
have been detected in the interstellar medium of the local 
galaxy [12] and some at high redshift [13]. Methylamine, 
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Figure 1. Variation of the potential energy of methanol as 
function of the relative rotation 7 of the OH-group with re- 
spect to the methyl group about the molecular axis. Shown 
are the J = 1, = 1 energies of the lowest torsion- 
vibrational levels. The splitting between the different sym- 
metry levels is due to tunneling through the potential bar- 
riers. The y4-symmetry species are split further due to the 
asymmetry of the molecule (if-splitting). 



another relatively small internal rotor molecule is compu- 
tationally more complex and will be treated in a separate 
paper [15]. 

This paper is organized as follows. In Section II, we 
will give a brief review of the theory of internal rotor 
molecules and outline how the torsional-rotational lev- 
els are numerically calculated using the belgi code [16]. 
In addition, we present approximate expressions for ob- 
taining the torsional energy splitting as function of the 
barrier height and the reduced moment of inertia, and 
compare this to the output of numerical calculations of 
a full Hamiltonian. In Section III, we will discuss how 
the molecular constants that appear in the torsional- 
rotational Hamiltonian scale with /z. These scaling re- 
lations are then used to determine the sensitivities of se- 
lected transitions in five different internal rotor molecules 
using BELGI. In Section IV the analytical expressions for 
the torsional energy splitting presented in Section II are 
used to construct a simple model for obtaining iiT^ from 
the three rotational constants [A, B and C) and three 
torsional constants (F, and p). This model provides 
an intuitive picture of the physics involved and makes it 
straightforward to estimate the sensitivity of other inter- 
nal rotor molecules. 



A. Hamiltonian 

The potential energy of an internal rotor molecule is 
a periodic function of the torsional angle 7 between the 
Csi, group and the remainder of the molecule, as shown 
in Fig. 1. Hence, it can be expanded in a Fourier series 
as 
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where typically Vg is about one hundred times smaller 
than V3, but provides information on the shape of the 
torsional potential. If only the first term of the expan- 
sion is taken into account, the torsional wave functions 
and energies follow from the solutions of the Mathieu 
equation [11]. 

To express the kinetic part of the Hamiltonian an axis 
system (a, &, c) must be defined. The symmetric part 
of the molecule is defined as the internal rotor or top 
whereas the remainder of the molecule is referred to as 
the frame, which for all molecules presented here has a 
plane of symmetry. The origin of the coordinate system 
coincides with the center of mass of the molecule. The a 
axis is chosen parallel to the symmetry axis of the top, 
and the b axis lies in the plane of symmetry. The c axis 
follows from the definition of a right-handed coordinate 
system. The inertia tensor then takes the form 



(3) 



with la = '^i 'nii{h'j -l-cf) being the moment of inertia 
of the molecule about the a axis. The subscript i runs 
over all atoms with mass rrii in the molecule. Expressions 
for If) and Ic, the moments of inertia around the h and 
c axis, respectively, can be found by cyclic permutation 
of the a, h and c labels. lab — '^iO'i^i is the product 
of inertia about the a and b axis. The kinetic energy can 
be expressed as [11] 
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lab^aUJb + 7;Ia2j^ + Ia2i^aj, (4) 



II. HINDERED INTERNAL ROTATION 

A review of hindered internal rotation can be found in 
the seminal paper by Lin and Swalen [11], while a recent 
review of various effective Hamiltonians, methods and 
codes dealing with asymmetric-top molecules containing 
one internal rotor with (or close to Csu) symmetry 
can be found in the paper by Kleiner [17]. In this sec- 
tion we will summarize those results that are relevant for 
obtaining the sensitivity coefficients. 



with LOa, i-Ob, and uJc the angular velocity components 
around the a, b, and c axis, respectively, and Ia2 the mo- 
ment of inertia of the top along its own symmetry axis. 
For a given vibrational state the zeroth order torsion- 
rotation Hamiltonian can be separated into a symmet- 
ric top part, an asymmetric top part and a torsional 
part [18, 19] 
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Table I. Some low-order symmetry-allowed torsion-rotation "Rho-Axis Method" (RAM) Hamiltonian terms for an asymmetric 
top containing a Csv internal rotor and a partial parameter list as used in the belgi code (taken from Kleiner [17]). The 
/Lt-dependence of the molecular constants is given in parenthesis. 







Torsional/Potential 






1 


1 — cos 37 




PaP-y 


1 — COS 67 




PaP% 


Rotational 






(^*°) 






(/^°) 






1 






V3/2 


F 


P 


V6/2 




k3 


p2 




(B + C)/2 


Fv 


Gy 


i„ 


Ny 


My 


ksj 


Pa 




A-(B + C)/2 


fes 


fe2 


fel 


K2 




ksK 


Pb - Pa 




{B~C)/2 


C2 


Cl 


C4 


cii 


C3 


C12 


PaPb+PbPa 




Dab 


dab 




Sab 


dabe 




SSab 



where 

His = 1{B + C) {P^ + P!) + API (5b) 

^^RA =\{B-C) {Pi - Pi) + Da, {PaPb + P^Pa) 

(5c) 

Hl^^ = F{p^+pPaf + V{^). (5d) 

Pa, Pb and Pc are the usual angular momentum operators 
along the a, h and c axis, respectively, and = —id/d^ 
is the angular momentum operator associated with the 
internal rotation of the top with respect to the frame. 
The coupling between the internal rotation and overall 
rotation in Eq. (5d) can be eliminated partly by trans- 
forming to a different axis system, the so-called "Rho- 
Axis System". In the resultant "Rho-Axis Method" 
(RAM), which is implemented in the BELGI code used 
hereafter, the torsional Hamiltonian operator contains 
only the +2FpPaP-y term. It is important to note here 
that two sign conventions exist in the literature for the 
torsion-rotation operator in Eq. (5d), i.e., F(py + pPa)^ 
and F{pj — pPa)^. If the latter convention is adopted the 
±K labeling of the E levels {vide infra) is reversed [20]. 
In this paper we adopt the convention with the "+" sign 
-i.e. F{pj+pPa)-^. The effective rotational and torsional 
constants are defined by 
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where lai is the moment of inertia attributed to the 
frame, defined by lai = Ia — Ia2- A dimensionless param- 
eter p is introduced by the axis transformation described 
earlier. For a symmetric top, p is simply defined as the 



ratio between the moment of inertia of the top divided by 
the moment of inertia of the molecule along the o-axis, 
i.e. p = Ia2/Ia- For asymmetric molecules p is a more 
complicated function of the various moments of inertia: 

p= '°;f^ - (11) 

^ah - iab 

The Hamiltonian of Eq. (5) can be expanded by adding 
additional distortion and interaction terms. Many of 
these higher-order constants multiply torsional operators 
by rotational operators and can be considered as effec- 
tive constants after the van Vleck transformations of the 
torsion-rotation Hamiltonian [17, 19, 21]. Some low-order 
symmetry-allowed torsion-rotation terms for an asym- 
metric top containing a C3,, internal rotor are listed in 
Table I. 

The overall Hamiltonian can now be written as 

H = H° + Hi^t + H^.i., (12) 

where Hc.d. corresponds to the centrifugal distortion 
Hamiltonian, and iJint contains higher-order torsional- 
rotation interaction terms. 



B. Eigenfunctions and Eigenvalues 

1. Torsion 

Hcrbst et al. [19] suggested to evaluate the Hamilto- 
nian of Eq. (12) in two steps. In the first step, the tor- 
sional Hamiltonian H^^^^ is diagonalized in a product ba- 
sis set composed of free rotor torsional eigenfunctions of 
p^ and eigenfunctions \K) = (27r)~^/^ exp(— i/Cx) of Pa, 
where x is the Euler angle: 

\Kka) = ^=\K)exp{i[3k + a]-f), (13) 
V 27r 

where a can take the values —1, 0, or +1 and k can be 
any integer. The eigenvalues of p-y are 3k + a as required 
by the periodicity of the potential. Due to the symmetry 
of the torsional Hamiltonian, basis bmctions of different 
a do not mix. Moreover, basis functions of different K do 
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not mix either, as Pa and H^^^^ commute. The resulting 
Hamiltonian matrix for each value of a and K is infinite 
in size, but it was found that truncating it to a 21 x 21 
matrix (— 10 < fc < 10) is sufficient to obtain experimen- 
tal accuracy for the molecules under study here [19]. The 
torsional eigenfunctions can be written as 
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exp(i [3fc + (T]7), (14) 



where vt is the torsional vibration quantum number and 
^3fe+CT ^''^ expansion coefficients. States with a — are 
labeled as A, and states with a — +1 and a — —1 are la- 
beled as El and i?2, respectively. For A torsional states, 
■izK levels are degenerate, whereas for E states a degener- 
acy exists between El, K and E2, —K levels. Although 
the torsional El and E2 state have different labels, tran- 
sitions between these two states are allowed. It was there- 
fore suggested by Lees [22] to refer to El and E2 levels 
as E levels where the sign of K distinguishes the two 
symmetries. 

The eigenvalues of the torsional Hamiltonian for 
methanol, acetaldehyde and acetic acid are depicted in 
Fig. 2 as a function of K for the ground torsional state 
[vt = 0) of those molecules. The solid circles, open cir- 
cles and open triangles represent values numerically cal- 
culated using the belgi code [16, 23-25] for the A, El 
and E2 torsional states, respectively. It is seen that the 
torsional energies are periodic functions of K with a pe- 
riod that is proportional to p^^ {p = 0.81, 0.33 and 0.07 
for methanol, acetaldehyde and acetic acid, respectively). 

In order to obtain an analytical model for estimating 
the sensitivities of transitions in internal rotor molecules, 
to be discussed in Section IV, we will now derive approx- 
imate solutions to the torsional Hamiltonian. It is clear 
from Eq. (14) that substituting a with cr -t- 3 results in 
identical eigenvalues, consequently, the eigenvalues may 
be regarded as periodic functions which can be expanded 
in a Fourier series as [11] 
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Figure 2. Torsional energies obtained with belgi for A 
(solid circles), El (open circles), and E2 (open triangles) lev- 
els as function of K for methanol (CH3OH), acetaldehyde 
(CH3COH), and acetic acid (CH3COOH) for ut = 0. The 
solid curves are fits to Eq. (15) for A, El, and E2 states. 
Note that only integer values of K have physical meaning. 



are plotted as the open diamonds in Fig. 3. By diago- 
nalizing H^^j.^ for several values of s and p > 0, while all 
other constants are set to zero, and fitting the torsional 
energies according to Eq. (15), ai coefficients were ob- 
tained for each value of s. These generic coefficients are 
plotted as the solid diamonds in Fig. 3. According to 
Lin and Swalen [11] the ai coefficients are given by the 
following equation; 



where ao,ai and higher order terms are (dimensionless) 
expansion coefficients. It can be shown [11] that these 
coefficients are functions of the reduced barrier height s, 
with 



9F ' 



(16) 



and that in the moderate to high-barrier limit, the series 
converges quickly. The solid curves shown in Fig. 2 are 
obtained by fitting the first two terms of Eq. (15) for A 
(cr = 0), El (cr = +1), and E2 (cr = -1) states. The re- 
sulting coefficients for methanol, acetaldehyde and acetic 
acid, as well as those for acetamide and methyl formate 
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(17) 



The solid line shown in Fig. 3 is obtained by fit- 
ting Eq. (17) to the generic ai values, using Ai = 
-5.296, Si = 1.111 and Ci = 2.120. Note that these 
fit parameters deviate from those given in Table IV of 
Ref. [11], but the curves agree with the curves shown 
in Fig 7 of the same paper. Small differences between 
the curves and the ai coefficients obtained for the differ- 
ent molecules can be attributed to higher order torsional 
terms which were not taken into account to obtain the 
fits. 
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Figure 3. The a\ expansion coefficients of Eq. (15) as a func- 
tion of tlie reduced barrier heiglit s for tlie ground torsional 
state ut = Q and first excited torsional state i^t = 1. Solid 
diamonds and squares represent ai coefficients in Ut = G and 
vt = 1, respectively, determined by fitting the eigenvalues of 
Eq. (5d) obtained with belgi for different values of a accord- 
ing to the expansion of Eq. (15). The solid and dashed curves 
are fits according to Eq. (17) and Eq. (29), respectively. Open 
diamonds and squares are ai expansion coefficients for several 
molecules taking into account higher-order torsional parame- 
ters. The inset shows an enlargement of the i^t = curve near 
the values for the six isotopologues of methanol. 



radiative transitions beween A and E species do not oc- 
cur. The arrangement of energy levels within a sym- 
metry state is quite similar to the structure of the K 
ladders in a (prolate) symmetric top, however in the 
case of internal rotation each K ladder attains an ad- 
ditional offset, Eto^s{K), induced by the tunneling split- 
ting. The -I- or - component of the A state refers to the 
I J, K, vu 0) ± I J, -K, uu 0) and | J, K, ut, 0) T I J, -K, 0) 
linear combinations of basis functions for K even and K 
odd respectively. The overall parity of the levels, is given 
by ±(—1)'^+''* [19]. When transitions between levels of 
different torsion-vibration states are ignored, the selec- 
tion rules for allowed transitions are [19] 



A levels: A J = |Aif|= (X 7^ 0), 1 

±o± |AJ|=1 |Aif|=0,l 

E levels: A J = |Ais:|= 1 

|AJ|==1 |Aif|==0,l 

where K is only a good quantum number in the limit 
of a symmetric top. As a consequence, transitions with 
\^K\ > 1 are allowed in asymmetric top molecules. 



III. SCALING AND SENSITIVITY 
COEFFICIENTS 



2. Rotation 

At this point the torsional Hamiltonian is diagonalized 
and the first step of the approach by Herbst et al. is 
complete. The second step of the approach consists of 
evaluating the remainder of the full Hamiltonian, i.e., 
overall rotation and coupling terms, in the basis set [19] 

\JKvt(j) ^\JK)\Kvta) , (18) 

with \ JK) the symmetric top rotational eigcnfunctions. 
The eigenvalues of the Hamiltonian of Eq. (5b) are 

{JKuta\H^,^,JJKvtCj) = 

]^{B + C)J{J+1)+[a-?-^^K\ (19) 

Note that, in the case of an asymmetric top, torsional A 
levels are split in ± components ioi K > Q due to the 
asymmetry of the molecule. 

C. Level Schemes and Selection Rules 

In Fig. 4 the lowest energy levels of i^CHj^OH are 
shown for the A and E species. A and E symme- 
try species can be considered as two different molecular 
species in the same sense as para- and ortho-hydrogen: 



The sensitivity to a variation of /i of a transition be- 
tween states [i/f, J", K", Ts") and J', K' . Ts) is given 
by 

K,{u':,J",K",Ts" -^u[,.r,K',Ts') = 
fiidE/dfj.)^,, j„ j^„^^^„ - f^{dE/dn)^, j,j^, r^^, 

E{iy'/, J", K", Ts") - E{v[, J', K' , Ts') ' ^ ' 

Thus, in order to calculate the coefficients, the energy 
of each level and its dependence on /i has to be obtained. 
This translates into knowing the values of the molecular 
constants that go into belgi and how these constants 
scale with /i. 

We will first examine the scaling relations for the lowest 
order constants. We implicitly assume that the neutron- 
electron mass ratio follows the same behavior as the 
proton-electron mass ratio, and that no effects depend- 
ing on quark structure persist [26]. As a consequence, the 
atomic masses and hence the moments of inertia are di- 
rectly proportional to /i. The rotational, centrifugal and 
torsional constants A, B, C, Dab, and F and the factor p 
are explicit functions of la, h, Ic, and lab, their /i depen- 
dence is obtained from Eq. (6)~(11). Within the Born- 
Oppenheimer approximation the torsional potential V3 is 
independent of the mass of the nuclei and hence of /x. 

The scaling relations for higher order constants are 
derived from multiple combinations of lower-order tor- 
sional and rotational operators. Some symmetry-allowed 
torsion-rotation terms are listed in Table I. Let us, for 
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Figure 4. Energy of the lowest rotational levels in the torsion- vibrational groundstate {i^t = 0) of methanol ( CHg OH) [23]. 
The levels are denoted by J (indicated on the left side of each level), K. For the A levels the so-called parity quantum number 
{+/—) is also used. The panel on the left contains the levels for the A state, whereas the panel on the right contains the levels 
for the E state. High sensitivities are expected for transitions that connect near degenerate levels with different K. 



example, inspect the constant M^, which can be consid- 
ered as a product of the torsional operator with the 
rotational operator P^. As the term with scales as 
and the term with scales as fJ,^^, we expect M^, 
to scale as The scaling relations of constants asso- 

ciated with other operators follow in a similar manner. 
The supplementary material of Ref. [10] lists the scaling 
relations for all constants used for methanol. 

In order to determine the sensitivity coefficients, we 
have written a computer code that generates the molec- 
ular constants as a function of fi using the discussed scal- 
ing relations, calls BELGI with these constants as input, 
and subsequently stores the computed level energies. 

As an example, the upper panel of Fig. 5 shows the en- 
ergies of the 22 E, 2i E,3^i E and 2o E levels in methanol 
as a function of A^/iJ,. The sensitivity coefficient of the 
2o — 3_i E and 2i — ^ 22 transitions can be obtained 
by dividing the difference in slope by the difference in 
energy, cf. Eq. (20). The sensitivity coefficients of these 
transitions as a function of the relative variation of the 
proton-electron mass ratio, A/i//i, are shown in the lower 
panel of Fig. 5. As expected from Eq. (20), the sensitivity 
is strongly enhanced when the energy difference between 
the levels becomes small. 

The behavior of the sensitivity coefficient close to the 



resonance depends on the quantum numbers of the lev- 
els involved. For instance, in methanol, K is not a good 
quantum number due to the asymmetry of the molecule, 
and levels with equal J mix. As a consequence, the 22 E 
and 2i E levels shown in Fig. 5 display an avoided cross- 
ing, and the sensitivity coefficient for the 2i — > 22 E 
transition is zero at the resonance. In contrast, J is a 
good quantum number and levels of different J will not 
mix. As a consequence, the sensitivity coefficient for the 
2o — > 3-1 -E transition becomes infinite at the resonance. 
In practice, we are only interested in the value of at 
A/i//i < 10~^ and the effects of avoided crossings is rele- 
vant only if the levels cross extremely close to Afi/ fi — 0. 
In our study we have only come across one transition 
that has a significantly reduced sensitivity coefficient as 
a result of mixing of the energy levels involved, namely 
the 8o -J> 8-iE transition in ^^CDJ^OD at 4.2 GHz with 
Kf^ — 0.7. Note that the sign of the sensitivity coef- 
ficients at A/i//i = is positive if the levels cross at 
A/x//x < 0, and negative if the levels cross at A/i//i > 0. 

Using the recipe described above, the sensitivity co- 
efficient of any desired transition in a molecule contain- 
ing a C^v symmetry group can be calculated. We have 
calculated the sensitivity coefficients of many (> 1000) 
transitions in methanol, acetaldehyde, acetamide, methyl 
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Table 11. Selected transitions [27] and Kfj, coefficients for different molecules, calculated with belgi (fourth column) using the 
constants from Refs. [23-25, 28, 29] and the approximate model (fifth column) discussed in the text. For methanol the most 
sensitive lines are shown, whereas for the other molecules the eight lowest transitions that have been observed in the interstellar 
medium as listed in the review by Lovas [27] are given. The error in the last digit(s) is quoted within brackets. Molecules 
marked with an asterisk (*) are labeled according to the sign convention opposite to the one as used in the text to be consistent 
with literature - i.e. F{p^ — pPa)^ instead of F{j)~f + pPaY ■ As a consequence the sign of the ±K labeling of the E transitions 
is negated for these molecules. 
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lo - 


► OoE 


19 262.140(4) 


— 1.00(5) 


— 1.00 




lo - 


► O0A+ 


19 265.137(1) 


— 1.00(5) 


— 1.00 




2o - 


► loE 


38 506.035(3) 


-1.00(5) 


-1.00 




2o - 


> I0A+ 




— 1.00(5) 


— 1.00 




1-1 


loE 


47 746.980(5) 


— 1.03(5) 


— 0.93 




li - 


► lo^+ 


47 820.620(4) 


— 1.02(5) 


— 1.03 


Acetamide* 


22 — 


- 2iA'^ 


9 254.418(4) 


— 1.04(5) 


— 1.12 




In — 
-^0 


- liJ5 


1 "^fis 7n^?A^ 

J-O 000. ( UO^^^ ) 


— 1 57(8) 


— 1 34 






, 42^T 


14 ZlU.o4y(4J 


— 1.03(5) 


— 1.12 




33 _ 


> 82^^ 


14 441.705(4) 


— 1.05(5) 


— 1.12 




2o — 


> 2i_E 


15 115.748(4) 


— 1.43(7) 


— 1.30 




22 - 


► liB 


zz uyo.oz i y^j 


— 74(5) 


Q 




33 - 


► 32A± 


ZZ i'Dy.DOC>[^4 J 


— 1.03(5) 


— 1.08 




42 — 


► 03-n/ 




— 1.04(5) 


— 1.11 


Methyl Formate* 


ii - 


► 1-1^; 


1610.900(2) 


-1.00(5) 


-0.70 




2-1 


-*■ l-iB 


22 827.741(8) 


-1.00(5) 


-1.00 




2i - 


>iiA+ 


22 828.134(8) 


-1.00(5) 


-1.00 




2o - 


> lo-E 


24 296.491(8) 


-1.00(5) 


-1.00 




2o - 


,UA+ 


24 298.481(8) 


-1.00(5) 


-1.00 




2i - 


<■ li_E 


26 044.796(8) 


-1.00(5) 


-1.00 




2i - 


> liA- 


26 048.534(8) 


-1.00(5) 


-1.00 




4o — 


- 3oE 


47534.069(16) 


-1.00(5) 


-1.00 


Acetic Acid* 


8-1 


-> 7-iE 


90 203.444(20) 


-1.00(5) 


-1.00 




8-1 


7oE 


90 203.444(20) 


-1.00(5) 


-0.97 




80 - 


► 7-iE 


90 203.444(20) 


-1.00(5) 


-1.03 




80 - 


7oE 


90 203.444(20) 


-1.00(5) 


-1.00 




80 - 


' 7oA+ 


90 246.250(50) 


-1.00(5) 


-1.00 




80 - 


► 7iA+ 


90 246.250(50) 


-1.00(5) 


-1.00 




81 - 


' 7oA+ 


90 246.250(50) 


-1.00(5) 


-1.00 




81 - 


► 7iA+ 


90 246.250(50) 


-1.00(5) 


-1.00 



formate and acetic acid, using the constants listed in 

Refs. [23-25, 28, 29], respectively. In Table II sensitiv- 
ity coefficients of selected transitions in tlie vibrational 
ground state (fj = 0) of these molecules are listed. For 
methanol, the most sensitive transitions involving lev- 
els with J < 10 are listed. For the other molecules no 
large sensitivities were found for the vibrational ground 
state and the eight lowest transition frequencies that 
have been observed in the interstellar medium have been 
listed [27]. Prom the table, we see in particular that 



transitions in methanol are much more sensitive to a 

variation of /i than the transitions listed for the other 
molecules. Except for the l_i — >■ li E transition in ac- 
etaldehyde with = —3.7, all transitions in acetalde- 
hyde, acetamide, methyl formate and acetic acid have 
-1.57 <Kf,< -0.74. 

The error in the last digit of the i^^-coefficients is 
quoted in brackets, and is conservatively taken to be 5% 
if \Kij,\ > 1 or 0.05 if jfC^j < 1. The error in the sen- 
sitivity coefficients has 3 sources: (i) errors due to the 
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Figure 5. Energies of selected rotational levels in methanol 
(^^CHg^^OH) as a function of the fractional change in the 
proton-to-electron mass ratio. The insets in the upper panel 
show magnifications of the regions where the levels cross. 
It can be seen that levels with same J give rise to avoided 
crossings, whereas levels with different J do not. The lower 
panel depicts the sensitivity coefficients of the 2i — > 22 -B and 
2o — > 3-1 E transitions as a function of Afi/ fi. For transitions 
between levels with same J, A'p goes to zero at the (avoided) 
crossing, whereas transitions between levels with different J 
Kfj, diverges to infinity at the crossing. 



uncertainty in the determination of the molecular con- 
stants. As the simulations reproduce almost all transi- 
tions <100 kHz, this error is negligible small, (ii) Errors 
due to inexactness of the scaling relations of higher order 
constants. Many of the higher-order constants are prod- 
ucts of torsional and rotational operators and may also 
be fairly correlated. Therefore, the exact relationships 
between the higher order parameters and the moments 
of inertia (and masses) are not obvious. For methanol 
we have investigated the influence of the higher order 
paramaters by comparing coefficients calculated by 
scaling only the first 7 and the first 31 constants to Kf^ 
coefficients obtained by scaling all 119 constants. Sen- 
sitivity coefficents were found to typically agree within 
5% or 0.5% if 7 or 31 constants were scaled, respec- 
tively. The effect of the higher order terms is therefore 
expected to be small. This has been confirmed by an 
independent study by Levshakov et al. [30]. (iii) Errors 
due to neglecting the fi dependence of the torsional po- 
tential. Within the Born-Oppenheimer approximation 
the torsional potential, V3, is independent of the mass 
of the nuclei and hence of /i. It is known however that 
V3 does vary between isotopologues. For instance for 



for ^^CDg^^OD V3 w 362 cm- 1. A reliable model for 
this variation is not available. As a check, we have as- 
sumed to be a linear function of I^^d — Iaila2/Ia] 
Fs^ysC'CHgi^OH) - 19.4(/,ed - /rcdC'CHgi^OH)). As 
/rod is directly proportional to /i, this introduces a de- 
pendence in the potential. We found that the coef- 
ficients for '^^CHg^^OH calculated by including the linear 
scaling for V3 are typically 3% smaller than those ob- 
tained when the potential is assumed to be independent 

of /i. 



IV. "TOY" MODEL 

Although the numerical calculations described in the 
previous section yield the sensitivity coefficient for any 
desired transition, they provide limited insight. In this 
section we will devise a simple model which provides an 
intuitive picture of the physics involved and aids in the 
identification of other internal rotor molecules that pos- 
sibly exhibit large sensitivity coefficients. In this model, 
we neglect coupling between vibrational, rotational and 
torsional motion. In this case, the n dependence of the 
energy of a certain state ji't, J, Ts) can be written as 



^/^ / ut,.J,K,Ts 



'vib 



where K'','^ 



^rot 



and 



(21) 

are the sensitivities to a pos- 
sible variation of the proton-to-electron mass ratio of a 
vibrational, rotational, and torsional transition, respec- 
tively. We neglect vibrational excitation and use Eq. (21) 
to rewrite Eq. (20) as 



- Ai?t,nrs 



(22) 



with ASrot and AE'tors being the difference in rota- 
tional and torsional energy between the two energy lev- 
els involved, respectively. From this equation it immedi- 
ately follows that diverges for Ai^i-ot — — AE'tors and 



This implies that the highest sensitivi- 



^CH 



OH the torsional potential V3 



373 cm , and 



^rot / j^t 

ties are expected for transitions that convert overall ro- 
tation into internal rotation or vice versa. Furthermore, 
Eq. (22) indicates that the coefficients are propor- 
tional to the amount of energy that is cancelled. We will 
now derive approximate analytical expressions for AE'^ot j 
A£;tors, a:™* and AT*"''" using the results of Sec. II. 

The rotational energy is given by Eq. (19), from which 
it is straightforward to calculate the energy differences for 
different transitions. The results are listed in Table III. 
Note that these expressions are valid only for a nearly 
symmetric molecule. 

AT™* follows from the fj, dependence of the rotational 
constants. From Eq. (6)-(8) A, B, and C are inversely 
proportional to the moments of inertia that are propor- 
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Table III. Rotational and torsional energy differences and associated sensitivity coefficients. Expressions for differences in 
rotational energy are only valid in the limit of a symmetric top molecule. 



Rotation AJ = AK = ±1 [A - ^{B + C)] {1 ± 2K) -1 

AJ = ±1 AK = ±(B + C) [J + i(l ± 1)] -1 

AJ = ±1 AK = ±1 ±(B + C) [J+ i(l± 1)] + [A- i(B + C)] (1±2X) -1 

AJ = ±1 AK = Ti ±(B + C) [J+ i(l± 1)] + [A- i(B + C)] (1=F2K) -1 

Torsion = 0" ±V3Fai sin [pK ^ ^]) (Si - 1) - i Ci \/s 

Aif = ±1 =F2Fai sin (f p) sin {p (a: ± i) - (t}] (Bi-l)-iCiyi 



Note tliat the transitions involving different overall symmetry species of the torsional levels [A E) are forbidden, thus transitions 
with /AK = 0, AJ = do not occur 




Figure 6. Torsional energy splitting between A and E J = 
0,K = (solid line) and El and E2 J = 1, \K\ = f levels 
(dashed line) as function of s from Eq. (23) and (24). For 
the curve representing the splitting between the El and E2 
states, p is fixed at p = 0.8. Also shown are the torsional 
splittings for the different isotopologues of methanol. 



As transitions between different symmetry states {A o 
E) are not allowed, we are interested in transitions within 
the same torsional states that differ in K. The difference 
in torsional energy for such transition with AK = ±1 
can be derived from Eq. (15) to be 

Actors = T2Fai sin (f p) sin (f {p (A- ± i) + a}). 

(25) 

Due to the oi-coefficient appearing in Eq. (23)-(25) 
all torsional splittings have the same dependence on s, 
and hence on p. The sensitivity of a torsional transition, 
K^J"^^, can be obtained from Eq. (1) by writing 



^tois _ 



/ aA^tors \ (ds\ fl 
\ ds J\dfiJ A£;tors 



ds 



AEu 



1, 



(26) 



where we have used the fact that F scales with 
The —1 appearing in the second line is introduced by 
the substitution of Actors with AEto^s/F- By inserting 
Eq. (25) in Eq. (26) we obtain 



tional to /i. Consequently, in first order approximation 

^lot ^ _i 

The torsional energy is given by Eq. (15), from which 
the splitting between the A (cr = 0) and E [a = ±1) 
states can be calculated 

AEtt"" = Fa, [cos {^-f{pK±l)]- cos {^-fpK)] 

^ ±^/iFai sin [pK ±\]). (23) 

Analogously, the splitting between El [a = +1) and E2 
(tJ = —1) levels is given by 

A<ir = ^Fa, sin {^-fpK). (24) 

Equations (23) and (24) are plotted in Fig. 6 as func- 
tion of s for i^T = and \K\ = 1, respectively, using 
the Ai, Bi and Ci parameters obtained in Sec. II. Also 
shown are data points for the various isotopologues of 
methanol derived from experiments. Within the Born- 
Oppenheimer approximation s is proportional to I^^d — 

\n'F-\ 



^tors ^ (^^ _ ^) _ 1^^^ 

~ 0.111 - 1.060^/s, (27) 

where we have used the dimensionless fit values for Bi 
and Ci obtained in Sec. II. Hence, within our approxi- 
mations, is only a function of s. For ^^CHg^^OH, 
with s = 6.01, this results in K^°^^ = —2.5. In agreement 
with the value found in Jansen et al. [10]. 

With the help of Eq. (22) and the expressions for 
AE.ot, AE^tors, if™' and AT*"" as listed in Table III, we 
can now determine which transitions are likely to have an 
enhanced sensitivity and estimate the coefficients of 
these transitions. From Eq. (22), we saw that the high- 
est sensitivities are expected when AE-^ot — ~ Actors ■ In 
Fig. 7 the difference in rotational energy, Ai?i-ot, (solid 
curves) and negated torsional energy, —AEy-ot (dashed 
curves) are shown for \J,K) \ J + 1, K — 1) yl+ transi- 
tions in methanol. The highest sensitivities are expected 
when the lines representing the difference in rotational 
energy AEj-^t and the negated difference in torsional en- 
ergy — AE'tors cross. For this to happen at low J and K 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 

J 



Figure 7. Difference in rotational energy, AiJrot, (solid 
lines) and negated difference in torsional energy, — AiStors, 
(dashed lines) between levels connected by a | J, K) — > 
|J + 1,K — 1) type transition in methanol. The curves 
are obtained by the expressions in Table III and the molec- 
ular constants of methanol [23]. The highest sensitivities are 
expected when torsional energy is converted into rotational 
energy or vice versa - i.e. when the two curves cross. The 
shaded area represents the maximum torsional energy that 
can be attained by the molecule. Note that the amount of 
energy that is cancelled is proportional to the Ky, coefficient. 



it requires that the rotational constants A, B and C are 
of the same order as the difference in torsional energy. If 
the rotational constants are much smaller or larger than 
Ai?t,ors, the crossings will only occur for high J and/or K 
quantum numbers. With the help of Fig. 7, it is straight- 
forward to select transitions that are likely to have a large 
coefficient. For instance, the solid line representing 
the difference in rotational energy for K = 1 (solid cir- 
cles), crosses the dashed line representing the negated 
difference in torsional energy near J — b. As the lines 
cross near the maximum of torsional energy that may be 
attained, represented by the border of the gray shaded 
area, we may expect a large A'p for the 5i — > 6o tran- 
sition. Indeed, this transition has = —42. 

The sensitivity coefficients K^^ of the transitions can 
be estimated using Eq. (22). Unfortunately, we found 
that the agreement between the coefficients obtained 
from this simple model and the values found from the 
full calculation was unsatisfactory, mainly as a result of 
neglecting the assymmetry of the molecules. Hence we 
chose to use the experimental energy difference between 
the levels, hv, rather than (Ai^rot + A_Etors)- In this case 
Eq. (22) can be written as 



= if™'+Ai?tor. {Kl°^^ - Kl°') /hv 

= -1 T iFAis^^e'^^-^' {Bi - x 
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s = 4rj9F 

Figure 8. The product f{s) ■ sin (f p) (see text) in the ground 
torsional state as a function of the effective barrier height, s, 
using Eq. (28) with p as indicated in the figure. Also shown 
are data points for each molecule investigated in this paper, 
values of s and p can be found in Table IV. 

sin (f p) sin {p {K±^)+ a])/hv 
= -l + F- fis) ■ sin (f p) • gip, K, a)/hy, (28) 

The calculated coefficients from this model for se- 
lected transitions are listed in Table II and are seen to 
agree rather well with the numerical calculations. 

We are now ready for a qualitative discussion of 
the sensitivity coefficients obtained for the different 
molecules. In the last line of Eq. (28), we have sepa- 
rated the expression for the sensitivity coefficient in four 
parts. The molecular constant F, a function /(s) that 
depends only on s, a function that depends only on p 
and a function, g{p^ K, a) that depends on the rotational 
quantum number, K, on the torsional symmetry, cr, and 
on p. This last function takes on a value between -1 and 
1 for the different cr, K levels. Although this function 
determines the sensitivity of a specific level, it is not im- 
portant for comparing different molecules. The product 
/(s) -sin (§p) can be used as a means to compare the sen- 
sitivity for different molecules. In Fig. 8 this product is 
plotted as a function of s with p as indicated in the figure. 
The curves can be regarded as the maximum sensitivity 
one may hope to find in a molecule with a certain F and 
transition energy hv, if g{p, K, a) is set to ±1. The max- 
imum sensitivity peaks at s = 4 and p — 1. Recall that 
p is defined as the moment of inertia of the top over the 
moment of inertia of the whole molecule (p ~ Ia2/Ia)j 
and cannot be greater than unity. 

In Table IV the structure and lowest order molecular 
constants of the molecules investigated in this paper are 
listed as well as the results from our analytical model. 
The last column lists the generic maximum sensitivity 
K^°'^ that may be expected for a hypothetical transition 
with a frequency of 1 GHz. For methanol, this number 
is 10-100 times larger than for the other investigated 
molecules, following from the fact that methanol has a 
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Table IV. Structure and some lower-order constants of the molecules investigated in this paper. The A — E K = torsional 
splitting from Eq. (23) is listed for the vibrational ground state i/* = of these molecules. Note that the magnitude of the 
torsional splitting only depends on the reduced barrier height s = 4V3/9F. The magnitude of Kj^'^^ for each molecule follows 
from Eq. (26). The K^'^'^ values given in the last column of the table are hypothetical sensitivities that may be expected for 
transitions of 1 GHz and with g{p,K,a) = 1. Molecular constants are taken from Refs. [23-25, 28, 29]. 





Structure 


Isotopo 


loguc 


V3 


F 


s 


P 






ifs- + 1 










(cm-i) 


(cm-i) 






(cm-i) 








H H 




OH 


373.555 


27.647 


6.01 


0.810 


9.07 


-2.5 


±398 








OH 


373.777 


27.642 


6.01 


0.810 


9.06 


-2.5 


±397 


Methanol 


/ 

c 




OH 


374.067 


27.428 


6.06 


0.809 


8.88 


-2.5 


±392 




/ 




OH 


370.055 


24.994 


6.58 


0.895 


7.12 


-2.6 


±363 




H 




OD 


366.340 


17.428 


9.34 


0.699 


2.58 


-3.1 


±145 








OD 


362.122 


14.758 


10.91 


0.822 


1.54 


-3.4 


±110 


Acetaldehyde 


H H 
H^. / 




C^'^OH 


407.716 


7.600 


23.84 


0.332 


0.065 


-5.1 


±3.6 



\ 



o 

H 



Acetamide c ^^CH,^^C ^''O "NH, 25.044 5.617 1.98 0.068 4.85 -1.4 ±5.2 



o 

o 



\ 

c 

»4 / 

/ 



Methyl Formate o ^^CHg'^O ^^C ^"OH 370.924 5.490 30.03 0.084 0.017 -5.7 ±0.3 



H 



Acetic Acid c ^^CH„^^C ^''O ^*OH 170.174 5.622 13.45 0.072 0.34 -3.8 ±2.8 



H 



large F, an effective barrier close to the optimal value of 
4, and a relatively large p. 

Application of the Toy model to excited torsional 
states 

So far, wc have limited the discussion to transitions in 
the ground torsional state {yt = 0). In this paragraph we 
discuss the application of the model to excited torsional 
states. Transitions within excited torsional states arc 
both unlikely to be observed in the interstellar medium 
and of less interest for laboratory tests due to the added 
complexity of the experiment, hence, the following dis- 
cussion is intended for providing a more complete picture 
only. In excited torsional states the splitting between 
the different torsional symmetry levels becomes larger as 
tunneling through the torsional potential becomes more 
likely. This implies that the energy that can be cancelled 
also becomes larger, but, at the same time, the sensitivity 
coefficient of a pure torsional transition, Kj^'^^, becomes 
smaller. As we will see, the first effect is more important. 

In Fig. 3 the solid squares represent ai coefficients for 



the i^t = 1 state obtained by fitting the second eigenval- 
ues of the torsional Hamiltonian as function of K to the 
Fourier expansion of Eq. (15). The dashed line in Fig. 3 
represents a fit using an expression similar to Eq. (17) 
with an additional term, i.e., 

a5;*=i = Ais^i • (29) 

with Ai = 10.388, = 0.829, Ci = 1.108 and A = 

—0.058. The open squares, also shown in Fig. 3, are oi 
coefficients for the first excited torsional state of the five 
molecules investigated in this paper. The additional term 
modifies Eq. (28) only slightly. With the known coeffi- 
cients, Ai, Si, Ci and Di, we can again plot the parts 
that depend only on s and p, as was done for i/f = 
shown in Fig. 8. As compared to the ft = 0, the curves 
for the fi = 1 are broadened and the center of the peak 
is shifted from s = 4 to s = 7. Moreover, the generic sen- 
sitivity is 5 times larger at the peak. As a consequence, 
the generic sensitivity, i.e., the sensitivity for a transition 
with 1/ = I GHz and g(p, K, a) = 1, of molecules with un- 
favorable s and p in the ground torsional state, can be 
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large in the t^t = 1 state. For instance, acetaldchyde in 
the t't = 1 state has a generic sensitivity if^'^" + 1 = ±82 
as compared to K^'^ + 1 = ±3.6 in the torsional ground 
state. For methanol, the generic sensitivity will also in- 
crease, however as the torsional splitting in the Ut = 1 
state is much larger than the rotational splittings, res- 
onances are expected to occur only at high J and K 
values. Note that the sensitivities from transitions from 
the i^t = to the vt — ^ state will not be significantly en- 
hanced as compared to ordinary vibrational transitions, 
i.e. K. -0.5. 



V. CONCLUSION &; OUTLOOK 

In the present study we have demonstrated that transi- 
tions in internal rotor molecules that convert internal ro- 
tation energy into overall rotation energy of the molecule 
exhibit an enhanced sensitivity to a possible variation 
of the proton-to-electron mass ratio. We have calculated 
the coefficients for five, relatively small, internal rotor 
molecules that are of astrophysical relevance; methanol, 
acetaldehyde, acetamide, methyl formate and acetic acid. 
In addition to full calculations using advanced codes such 
as the BELGI program achieving spectroscopic accuracy 
on the level energies, we have developed an approximate 
model, dubbed as "toy" -model, in which the molecular 
structure is described by the six most relevant of the 
molecular parameters. Based on this model we produce 
insight in the question as to why certain molecules of C^y 
symmetry with hindered internal rotation are sensitive 
to /i-variation. In particular, molecules in the torsional 
ground state are expected to have large if^ coefficients if 
they have (i) an effective barrier height, s, around 4, (ii) 
a ratio between the moment of inertia of the top and the 
whole molecule, p, close to unity, and (iii) a large value 
for the molecular constant that relates to the internal ro- 
tation, F. If the torsional splittings are of the same order 
as the rotational constants, sensitive transitions will oc- 
cur between levels with low J and K quantum numbers. 

From the approximate "toy" model we learn that of the 
five molecules studied, methanol has by far the largest 
sensitivity, due to its favorable value of the effective bar- 
rier s, and the fact that in methanol p is near unity. 
Moreover methanol has a fairly large value of F. The 
other investigated molecules either have a too large bar- 
rier (acetaldehyde), have a too heavy frame attached to 
the methyl group and consequently a very small p (ac- 
etamide and acetic acid), or have both a high barrier and 
a small p (methyl formate) . Based on these criteria, other 
interesting molecules containing a Cat, symmetry group 
include mercaptan, CH3SH, [F — 15 cm~^, V3 = 439 
cm~^ and p = 0.65, resulting in K^'^ + 1 = ±67) and 
methylamine (which will be topic of a separate publi- 
cation [15]). Other interesting candidate molecules, al- 
though of C2 symmetry, are H2O2 (recently treated by 
Kozlov [31]) and H2S2. These molecules require a modi- 
fication in the definition of s and will have different Ai , 



Bi and Ci coefficients. 

The high sensitivities of internal rotor molecules, par- 
ticularly methanol, make them excellent target species 
to search for a variation of the proton to electron mass 
ratio over cosmological time scales. It is important to 
note that the sensitivity coefficients of the transitions 
in these molecules have both large positive and large 
negative values, i.e., if the proton-to-electron mass ratio 
varies, some transitions will shift to higher frequencies 
while others shift to lower frequencies. This makes it 
possible to perform a test of the variation of the proton- 
to-eleetron mass ratio using transitions in pertainig to 
a single species, thus avoiding the many systematic ef- 
fects that plague tests based on comparing transitions in 
different molecules. 

Currently, the most stringent bounds on a cosmological 
variation of p is set by observations of hydrogen molecules 
in high (z = 2 — 3) redshift objects [32, 33] and a com- 
parison between the ammonia tunneling frequencies and 
rotational transitions in anchor molecules at intermedi- 
ate redshift {z = 0.5 — 1) objects [34-36]. Methanol 
provides a system that should result in more stringent 
bounds for y^-variation. Recently, methanol, as well as 
methylamine and acetaldehyde have been observed in the 
gravitationally lensed system, PKS 1830-211 a,t z = 0.89 
[13, 14]. The Iq 2_i E transition in methanol reported 
by Muller et al. was calculated to have -ftT^ = —7.4. We 
have also calculated the ii'^ coefficients for the nine ob- 
served transitions in acetaldehyde and found that all lines 
have sensitivities of K^j^ ~ — 1. Sensitivity coefficients for 
methylamine, for which three lines were observed in the 
same survey, will be calculated in a separate paper [15]. 

The high sensitivity coefficients in methanol are also 
beneficial for probing variation of /i as a result of 
chameleon- like scalar fields. These fields predict a de- 
pendence of p on the local matter density. Note that 
the physical origin of these chameleon like theories are 
very different from theories describing temporal p vari- 
ation. Levshakov et al. [37] compared ammonia spectra 
taken at high (terrestrial) and low (interstellar) densities 
of baryonic matter, and observed a statistically signifi- 
cant variation of p. Recently, a preliminary tests using 
methanol was performed by Levshakov et al. [30]. This 
test obtained similar results as for ammonia but can be 
further improved if more accurate laboratory data be- 
comes available. 

Methanol is also a promising candidate for laboratory 
tests on a possible variation of p. In laboratory experi- 
ments rare isotopologues can be used, in contrast to cos- 
mological searches. Hence the most sensitive transitions 
deriving from the present calculations can be targeted; a 
combined measurement of the 2i li E and 3o ^ A-i E 
lines in CD3OH exhibit a sensitivity spread Aii'p that is 
more than 400 times larger than a pure rotational tran- 
sition [10]. Note that the large dipole moment and low 
mass of methanol make it possible to use advanced labo- 
ratory techniques for cooling and manipulating molecules 
by electric fields [38, 39]. 
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